Possible galactic sources of ultrahigh-energy cosmic rays and a strategy for their 

detection via gravitational lensing 
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If decays of superheavy relic particles in the galactic halo are responsible for ultrahigh-energy cosmic 
rays, these particles must be clustered to account for small-scale anisotropy in the AGASA data. 
We show that the masses of such clusters are large enough for them to gravitationally lens stars and 
galaxies in the background. We propose a general strategy that can be used to detect such clusters 
via gravitational lensing, or to rule out the hypothesis of decaying relic particles as the origin of 
highest- energy cosmic rays. 
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The origin of cosmic rays jl]H with energies beyond 
the Greisen-Zatsepin-Kuzmin (GZK) cutoff jl is un- 
known. One of the possible explanations invokes de- 
cays of metastable superheavy relic particles X with 
masses 10 13 GeV or higher and cosmologically long life- 
times Such superheavy particles could be pro- 
duced non-thermally at the end of inflation [pd|-p^[. 
Their extremely small decay width may be due to a con- 
servation of some topological charge |J. Particles with 
the requisite properties may also arise from string the- 
ory 

If these particles decay into hadrons and photons, the 
flux of ultrahigh-energy cosmic rays (UHECR) is dom- 
inated by those particles in the halo of our galaxy ||. 
This can explain the absence of the GZK cutoff. Even if 
the superheavy particles decay predominantly into neu- 
trinos flq ], cosmic rays with energies beyond the GZK 
cutoff may originate through Z-bursts ||l6f . In this letter 
we concentrate on the former possibility and assume that 
observed ultrahigh-energy events come mainly from 
the decays of relic particles in the Milky Way halo. 

The new data provide an opportunity to test this 
hypothesis through gravitational lensing. There is a 
strengthening evidence for directional clustering of events 
in the AGASA data ||[|. Tne latest analyses || show 
one triplet and six doublets of events, each originating 
from the same point in the sky, to ±1.3° accuracy. The 
probability of this clustering to occur by accident is less 
than 0.07% [|J. The only way to reconcile these data with 
the hypothesis of relic particle decays is to assume a non- 
uniform distribution of particles in the halo. If the relic 
particles form regions of increased density, such lumps 
may be responsible for the doublets in the UHECR data. 
To produce a doublet, a lump of particles must be of a 
certain size determined by the decay probability. Since 
the mass of the hypothetical particle is fixed by the en- 
ergy of UHECR, there is a prediction for the mass of each 



lump that can give rise to a doublet. In addition, the ce- 
lestial coordinates of the particle cluster are known to one 
degree accuracy. In this letter we propose a novel gravi- 
tational lensing technique that can be used to discover a 
cluster of relic particles, or to rule out such particles as 
the origin of ultrahigh-energy cosmic rays. 

Under the assumption that UHECR are caused by the 
relic particle decays, the data suggest that (in addition 
to a possible uniform distribution), of the order of ten 
clusters of X-particles exist in our galactic halo. N-body 
simulations of dark matter halos predict some inhomo- 
geneities that can be related to small-scale anisotropy of 
UHECR H but probably are not sufficient to produce 
larger clumps. However, additional interactions of the 
hypothetical particles can alter this picture dramatically. 
We assume that each doublet comes from a separate clus- 
ter of particles. Here we do not discuss the dynamics of 
clustering of the heavy particles. This issue will be ad- 
dressed in an upcoming publication. We note in passing 
that clumps of dark matter with masses 10 8 M Q may re- 
solve 0| the widely debated issue of cusps in the halo 
density profiles M. 



If the X-particlc lifetime is t x 
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cluster of N particles produces decays at a rate 
P= — ~ (10- 10 - HT 22 )iV yr" 1 . 



(1) 



The probability for the decay products to produce an 
air shower in a detector is P x (d/L) 2 , where d ~ 10 6 cm 
is the size of the detector and L ~ 10 23 cm is the distance 
to the cluster of relic particles. In order to have a doublet 
in a one-year data set, each cluster must have 
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particles. If X-particle has mass m x ~ 10 13 GeV §§, 
the mass of the cluster is 
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FIG. 1. A typical light curve for gravitational lensing of a 
single source as a function of time (a), and the corresponding 
plot of brightness variations for a sample of stars (b) as a func- 
tion of celestial coordinates, in arbitrary units. Photometric 
changes indicate the location of the lens. 
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The lifetime t x can be in the range from 10 10 yr (for 
the relic particles to survive until present) to 10 21 yr (for 
the total mass of the clusters not to exceed the mass of 
the galaxy). Correspondingly, the masses of clusters can 
range from one solar mass to 10 10 AfQ. 

There is a remarkable possibility to discover such invis- 
ible massive objects by what we will call a "lens-chasing" 
technique. Although the cluster can be entirely dark, it 
can be detected through gravitational lensing of stars and 
galaxies behind it. AG AS A data |Q,E| provide the celes- 
tial coordinates of the clusters with a precision of a few 
degrees. 

The Einstein radius of a cluster with mass M is 



R E = 0.14 pc 
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A lens of this kind has angular size 
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and passes the line of sight in time 
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We propose to use the small changes in brightness of 
stars and galaxies behind the lens on time scales of the 
order of one year. The AG AS A data specify the loca- 
tion of the cluster to ±1.3°. It is possible to scan over a 
large sample of remote sources in the patch specified by 
the UHECR data, recording the brightness of the back- 
ground stars and galaxies. The scan must be repeated af- 
ter a period of several months. Next, one should extract 
the changes in the absolute brightness of the background 
stars. A slowly moving lens with a single-star light curve 
shown in Fig. 1 (a) will produce a map of brightness dif- 
ferentials shown in Fig. 1(b). By using temporal changes 
in the brightness of stars from a large sample, one can 
locate a small lens within a large, 2.5°, patch of the sky. 

The number of background sources chosen for scanning 
and photometry determines the sensitivity of the pro- 
posed lens-chasing experiment. Let us consider a sam- 
ple of n 2 stars with an average angular separation of 
9b « 2.5° /n. A lens that passes near one of these stars 
at an angular distance 59 = 9^/2 < 9 E from the line 
of sight will magnify the source by a factor (29 e 1 9b) 2 as 
compared to its brightness in the absence of the lens. The 
change in the star's brightness over a one year period is 
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Assuming a better than 1% precision photometry, a 
sample of 10 7 background sources allows detection of 
clusters with mass 1O 7 M and higher. Smaller masses 
require a higher number of background sources. For 
comparison, MACHO project has monitored 11.9 million 
stars during 5.7 years of operation jL7|. Of course, lens 
chasing presents a very different challenge from that faced 
by MACHO. Unlike MACHO, which monitors bright 
nearby stars in Large Magellanic Cloud (LMC) on a con- 
tinuous basis, we want a relatively infrequent (once a 
year) accurate photometry of stars in the directions of 
UHECR doublets. 

Some of the clusters in the AGASA data lie in the 
supergalactic plane |^|j3|. The presence of many rela- 
tively close (and, hence, bright) background stars in these 
directions makes the corresponding clusters particularly 
appealing for lens chasing. 

One can, of course, refine this technique. If the lensing 
cluster is discovered after several initial crude scans, one 
can narrow down its coordinates and perform a more 
detailed monitoring of closely spaced sources around the 
location of the lens. 

We note in passing that future detectors can observe 
yet another signature of the same kind of sources. De- 
cays of superheavy particles and subsequent fragmenta- 
tion can produce excited hadrons. Their decays, in turn, 
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can produce simultaneous air showers separated by thou- 
sands of kilometers. The time delay is 5t = tj~ 2 AE/E, 
where t is the time of flight, 7 ~ 10 11 is the Lorentz fac- 
tor, and AE/E ~ 1. The difference in the arrival time 
St ~ 10~ n s, and the distance between air showers is of 
the order of 10 3 km. Future space-based detectors, such 
as EUSO and OWL, can observe such spatially separated 
events in coincidence. 

To summarize, a small-scale anisotropy in the AGASA 
data demands that, if the UHECR are due to decaying 
relic particles in the halo, these particles form clusters 
with coordinates specified by the cosmic ray events. The 
masses of such clusters can range from M Q to 10 M®. 
We have proposed a general strategy for detecting such 
clusters by their gravitational lensing of the background 
stars and galaxies. 

This work was supported by the NATO Collaborative 
Linkage Grant PST.CLG.976397. In addition, A.K. was 
supported by the US Department of Energy, grant DE- 
FG03-91ER40662, Task C. The work of V.A.K. was sup- 
ported in part by the grant 98-02- 1744a from Russian 
Foundation for Basic Research. V.A.K. thanks UCLA 
for hospitality during his visit, when this work was per- 
formed. 



[10] K. Hamaguchi, Y. Nomura and T. Yanagida, Phys. Rev. 

D59, 063507 (1999). 
[11] D.J. Chung, E.W. Kolb, and A. Riotto, Phys. Rev. Lett. 

81, 4048 (1998); Phys. Rev. D59, 023501 (1999). 
[12] V. Kuzmin and I. Tkachev, JETP Lett. 68, 271 (1998). 
[13] For review, see, e.g., V. A. Kuzmin and I. I. Tkachev, 

Phys. Rept. 320, 199 (1999). 
[14] K. Benakli, J. Ellis, and D. V. Nanopoulos, Phys. Rev. 

D59, 047301 (1999); Phys. Rev. D59, 123006 (1999); 

M. Birkel and S. Sarkar, Astropart. Phys. 9, 297 (1998). 
[15] G. Gelmini and A. Kusenko, Phys. Rev. Lett. 84, 1378 

(2000) ; J. L. Cro oks, J. O. Dunn, and P. H. Frampton, 

astro-p h/0002089| . 
[16] T. Weiler, Astropart. Phys. 11, 303 (1999); D. Fargion, 

B. Mele and A. Salis, Astrophys. J. 517, 725 (1999). 
[17] C. Alcock, R.A. Allsman, D.R. Alves et at [MACHO 

Collaboration], Astrophys. J. 542, 281 (2000). 



[1] M. Takeda, N. Hayashida, K. Honda et al., Phys. Rev. 
Lett. 81, 1163 (1998); M.A. Lawrence, R.J. Reid and 
AA. Watson, J. Phys. G G17, 733 (1991); D. J. Bird, S. 
C. Corbato, H. Y. Dai et al, Phys. Rev. Lett. 71, 3401 
(1993). 

[2] M. Teshima, talk presented at First International 
Workshop on Radio Detection of High-Energy Particles 
(RADHEP-2000), Los Angeles, CA, November 16-18, 
2000. 

[3] M. Takeda, N. Hayashida, K. Honda, et al., Astrophys. 
J. 5 22, 225 (1999); N Hayashida, K. Honda, N. Inoue et 



'., a,stro-ph/0008102 



[4] K. Greisen, Phys. Rev. Lett. 16, 748 (1966); G. T. Zat- 
sepin and V. A. Kuzmin, Pisma Zh. Eksp. Teor. Fiz. 4, 
114 (1966) [JETP Lett. 4, 78 (1966)]. 
[5] V. Berezinsky, M. Kachelriess, and A. Vilenkin, Phys. 

Rev. Lett. 79, 4302 (1997). 
[6] P. Blasi and R. S heth, Phys. Lett. B 486, 233 (2000) 
|a,stro-ph/0006316 | 



[7] A. El-Zant , I. Shlosman, and Y. Hoffman, astro- 
ph/0103386|. 



[8] D. N. Spergel and P. J. Steinhardt, Phys. Rev. Lett. 84, 
3760 (2000); R. Dav e, D. N. Spergel P. J. Steinhardt 
and B. D. Wandelt, fa,stro-ph/0006218| ; B. D. Wandelt, 
R. Dave, G. R. F axrar, P. C. McGu ire, D. N. Spergel and 
P. J. Steinhardt, |astro-ph/0006344| . 

[9] V. A. Kuzmin and V. A. Rubakov, Phys. Atom. Nucl. 
61, 1028 (1998) [Yad. Fiz. 61, 1122 (1998)]. 



3 



